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Abstract 

In this article we show that one can construct initial data for the Einstein equations 
which satisfy the vacuum constraints. This initial data is defined on a manifold with 
topology R 3 with a regular center and is asymptotically flat. Further, this initial data 
will contain an annular region which is foliated by two-surfaces of topology S 2 . These 
two-surfaces are future trapped in the language of Penrose. The Penrose singularity 
theorem guarantees that the vacuum spacetime which evolves from this initial data is 

future null incomplete. 
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1 Introduction 

This paper is the third in a series dealing with the existence of globally regular, asymp- 
totically flat, maximal initial data for General Relativity in vacuo which have trapped 
surfaces. Our approach consists of considering certain classes of sequences of initial data 
called critical sequences (CS's) which, for large values of the parameter, develop regions 
closely resembling the black hole region of the extended Schwarzschild (SS-) spacetime. 
We first recall (see [1,2]) the properties of the SS-geometry which form the model for our 
more general results. 

Consider the spacetime metric 



ds 2 = -(l-^)dt 2 +(l-^) dr 2 + r 2 d0 2 (1.1) 

where < 2m < f < oo, — oo < t < oo and dO 2 is the line element on the unit two-sphere. 

The spacetime (1.1) can be viewed as the right quadrant of the Kruskal manifold. In 

particular it can be smoothly extended to the black hole region (i.e. the upper quadrant of 

Kruskal) by replacing t by the Eddington-Finkelstein coordinate r = t+f+2m log \ ^ — 1| 

and allowing f to vary over < f < oo and — oo < r < oo. For f < 2m we can, if we 

wish, again represent ds 2 by Equ. (1.1). Now fix a constant C > with 

27m 4 . . 

C 2 < 1.2 

16 V ; 

and consider the function ip : R + — > R defined by 

2m C 2 . . 

y,(r) = l- — + -. (1.3) 

tp is negative for f = 3m/2, positive for f = 2m, goes to one as f — > oo and is monotoni- 
cally increasing for f > 3m/2. Thus, for some unique r c with 3m/2 < r c < 2m, (p(r c ) = 
and ip > for f > r c . Moreover r c is a simple zero of ip. Now define a function / (f > f c ) 
by 

f{r) = -C f _. (1.4) 

J rc s 2 {l-2m/s)^(s) 

This will be the height-function of a spherically symmetric maximal slice through the 
Schwarzschild solution of the kind discussed in [1,2]. If (1.4) is taken in the sense of the 
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Cauchy mean value at s = 2m, /(f) is denned for f c < f < 2m and 2m < r < oo. One 
can see that /(f) is everywhere positive and is logarithmatically divergent at f = 2m. 
One can also see that the function g given by 

g(r) = /(r) + 2m log |- l|+r (1.5) 

is C°° for all f e (f c , oo). Thus the height-function diverges logarithmically at the hori- 
zon in the standard Schwarzschild coordinates but is everywhere regular in the regular 
Eddington-Finkelstein coordinates. Now consider the subset (of the extended manifold) 
given by r = g{f). For each fixed C satisfying (1.2) this is a smooth submanifold along 
which, for f ^ 2m, one has t = /(f). This will lie in the upper half of the right-hand 
quadrant since we know that / > and continues through to the middle of the upper 
quadrant of the Kruskal diagram. We can express this as f = / _1 (t), t > 0. The latter 
function can be smoothly extended to negative values of t by writing f = / -1 (|£|). This 
results in an extended submanifold Mq which is a Cauchy slice of the Kruskal manifold 
and lies symmetrical with respect to the timelike cylinder given by t = 0. (Note that the 
set where t — consists of a time symmetric Cauchy slice plus a timelike cylinder which 
intersect orthogonally along the bifurcation sphere of the horizon.) The two asymptotic 
ends of Mq are translates under d/dt of one another by the amount At = 2 lim-^^ /(f). 
This surface Mq lies is the upper half of the Kruskal diagram. The 'mirror' surface in the 
lower half plane corresponds to choosing C < 0. 

One can now compute the first and second fundamental form of M C - One gets (f c < 
r < oo) 

g ab dx a dx b = -dr 2 + r 2 d0 2 (1.6) 

and 

p ab dx a dx b = C [ -^df 2 - \d0 2 ) , (1.7) 
\^p r i> r j 

where a, b run from 1 to 3. Note that our choice of the definition of extrinsic curvature 
starts from K^ v = +jC u ^ flu which agrees with Wald[3] but is the negative of the choice 
made in MTW[4]. It now follows that p = p a b9 ab = and that p ab . a = 0. There is a unique 
(up to a constant factor) spherically symmetric 2-tensor which is both divergence-free and 
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trace-free (TT-) on flat space [5]. TT-tensors are conformally covariant. The metric g a b 
above is conformally flat and the tensor p^ is the conformal transform of the flat-space 
spherically symmetric TT-tensor. In flat space the TT-tensor blows up at the origin, but 
since the slice through Schwarzschild has topology S 2 x R it is everywhere regular. 

The complete 3-manifold M c consists of two copies of the asymptotically flat manifold 
(1.6,7), glued together along the minimal surface f — r c where f vanishes. (The subsets 
f > 2m of these two copies are contained in the right and left quadrant of Kruskal, 
respectively.) Thus the M c 's are maximal spherically symmetric Cauchy surfaces. Since 
r c > 3m/2, they stay away from the singularity r = uniformly in C. It is well known 
that the spheres of constant t — to, r — tq for tq < 2m are future trapped surfaces 
(FTS's). For spheres lying in Mc for some C, we can express this fact in terms of g a b and 
Pab- The divergence of the null normals to a sub manifold E of a spacelike slice M is given 
in general by (note that our choice of sign for p ab plays a role here) 

e± = H TPa b n a n b ±p, (1.8) 

where + (— ) refers to the future (past) outer null normal, n a is the outer normal to 
E embedded in M and H the mean curvature of E w.r. to n a . "Outer", in the present 
context, is just conventional, say "right" . The surface E is called future-trapped iff 6 + < 
and ©_ > 0. In such a situation the Penrose singularity theorem [6] implies that any 
vacuum spacetime Cauchy evolving from M is future null incomplete. 

Suppose, contrary to the Schwarzschild situation, that E divides M into an outer, 
non-compact region (containing spatial infinity) and a compact inner region. In this 
case a is called an outer trapped surface (OTS) when merely 0+ < 0. (Note that the 
notion of OTS, as opposed to that of FTS, is not necessarily time-asymmetric. We could 
have surfaces E for which both + and ©_ are negative.) It is known (see [7]) that the 
OTS-condition is already sufficient for the Penrose theorem to hold. 

Let us remark at this point that the time-symmetric slices of SS have neither FTS's 
nor OTS's (at least no spherically symmetric ones), so the Penrose singularity theorem 
does not apply. But, due to the existence of two infinities (since the topology is S 2 x R), 
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the singularity theorems of Gannon [8] and Lee [9] can be used, again with the result 
that any Cauchy evolution is geodesically incomplete. 
Applying (1.8) to f = const in M c we find that 




for the right copy and 




(1.10) 



for the left copy. Since (p > C 2 /r 4 for f > 2m and (p < C 2 /f 4 for f < 2m, 6 + and 6_ arc 
both positive for f > 2m in the right copy. At f = 2m the quantity 6+ changes sign and 
©_ remains positive. As we cross f = 2m in the left copy, + remains negative and 0_ 
changes sign. In particular, all spheres of constant vq in Mc with ¥q < 2m are FTS's. 

All spherically symmetric, complete maximal spacelike slices of SS are either the 
Mc's and their translates under d/dt or the time symmetric slices t = const which can 
be viewed as limits of the M c 's as C goes to zero. On the latter slices one has + = 0_ 
positive on the right copy, zero on the bifurcation 2-sphere and negative on the left copy. 

n n 

In previous work [10,11] we have constructed sequences of solutions (g a b,p a i,) with 
p = of the Einstein vacuum constraints on R 3 which are asymptotically flat, complete 
and which, as n — > oo, develop annular regions A n foliated by 2-spheres for which 0+ is 
negative. Although these data are not necessarily time-symmetric their behaviour in A n 

n 

is dominated by g for large n and, in particular, 0_ will also be negative in A n . In fact, 
for large n, the geometry of A n is approximated by an annulus of the left (in keeping 
with the previous convention) copy of a t = slice of SS. For SS of course, the left and 
right copy are completely equivalent, but in the CS's we consider the interior of A n is 
topologically a ball rather than a second asymptotic end. Thus the (extended) Penrose 
theorem on OTS's applies to spacetimes evolving from these data, and one concludes 
that, for large n, they are null geodesically incomplete both in the future and in the past. 
In the present work we perform a similar construction, but where we are able to control 
the sign of the second term in Equ. (1.8) in such a way that, for large n, we obtain annular 
regions A n in our initial data sets which are modelled on the manifolds (Mc,g a b,Pab) of 
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SS described above in an annular region close to the minimal surface r = r c . In particular 
the A^s for large n are foliated by future trapped surfaces. Whether, as is suggested by 
the Schwarzschild case, there is inside of A n a region which is both outer future- and 
past trapped, is at present beyond our control. Ultimately we are interested in initial 
data with FTS's which result from Cauchy evolution of ones leaving neither OTS's nor 
FTS's, or, more ambitiously, which are asymptotically flat at past null infinity so that 
they, provided weak cosmic censorship is true, describe gravitational radiation collapsing 
to a black hole. To construct such data, or perhaps to show that the ones constructed 
here have this property, is an open problem. 

In our most recent work on this topic we considered initial data sets with non-zero 
extrinsic curvature. However, we assumed that the extrinsic curvature fell off rapidly at 
infinity and played no significant role in trapped surface formation. The key improvement 
in this current article is that we control the asymptotic behaviour of the extrinsic curva- 
ture. In particular, we assume that the extrinsic curvature near infinity is dominated by 
the spherically symmetric TT-tensor of Equ.(1.7). Thus we construct a family of initial 
data, in which, to leading order, the intrinsic geometry looks like (1.6) and the extrinsic 
curvature looks like (1.7) with some constant C. This means that, to leading order, the 
null expansions approximate Equ.(1.9). Further, as in our previous work, we can allow the 
ADM mass to become unboundedly large while controlling the error terms. This means 
that we can construct initial data for which the now approximate formulae Equ.(1.9) are 
accurate for f > 2m. This initial data set will contain an annular region around f = 2m 
of future trapped surfaces with a minimal surface in the annulus. 

2 Definitions and results from previous work 

Let M be a compact manifold diffeomorphic to S 3 . For g a smooth Riemannian metric 
on M, define the conformal Laplacian acting on functions by 




A 3 + ^%], 
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where A g = g ab D a D b , D a the covariant derivative and 1Z the scalar curvature of g. The 
elliptic operator L g , viewed as a densely defined operator on the Hilbert space L 2 (M,g) 
is essentially self-adjoint, with real eigenvalues bounded from below. Let Xi(g) be it's 
lowest eigenvalue. Suppose Xi(g) > 0. Let A be some arbitrary point in M. Then L g 
has a unique, positive Green function with source point A. More precisely, there exists 
G : M \ A ^ R, G > 0, satisfying 



/ 

Jm 



(L g f)GdV g =4nf\ A (2.2) 

I M 

for all / G C°°(M), or equivalently 

L g G=4n5\ A , (2.3) 

where dV is the Riemannian volume element and 5\ A the Dirac delta distribution with 
source point A. The singularity of G near A can be described as follows. Let f2 be an 
asymptotic distance function (ADF). This should mean that VL > outside A and 

n\ A = 0, D a Q\ A = 0, (D a D b n-2g ab )\ A = 0, (2.4) 

D a D b D c Q\ A = 0. (2.5) 
This, in a local coordinate neighbourhood centered at A, is equivalent to 

n = n + ^g ab ,Mx a x h x c + o°°(n*), (2.6) 

where fi = 9ab(0)x a x b and / = O°°(n k /2 ) mean / = 0(fij /2 ), Df = 0(Q { o~ 1)/2 ), DDf = 
0(Q { o~ 2)/2 ), a.s.o. In particular fl = O°°(fi ). A function of the form SI = tt + 0°°(ttl /2 ) 
satisfies (2.4) but not necessarily (2.5). Given an ADF, G has the property that 

G = n- 1 / 2 + ^ + o(n 1 / 2 ) (2.7) 

for a constant m and 

d(G-tt- l l 2 ) = 0(1). (2.8) 

The constant m has the interpretation of the ADM mass of the metric g ab = G A g ab on 
M \ A = R 3 . Let u be positive. Then 

L § $ = u- h L g <$> (2.9) 
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for all <£>, where 

g = u A g, $ = w _1 $. (2.10) 
Letting uj = $ and $ = G, so that <5 = 1, it follows from Equ. (2.3) that 

TZ[g]=0 on M = M\ A (2.11) 

where 5^ = G A g a b- By virtue of (2.4,5) and (2.7) is asymptotically flat near A. Thus, 
by the positive-mass theorem [12,13], we have that m > and m = implies that g is 
the flat metric on M = R 3 , or, equivalently, <? is conformal to the standard metric on S 3 . 
Next let p be smooth in M \ A, the singularity at A being restricted by 

p = 0(^o 3 ), dp = O(fi 7/2 ) (2.12) 

and consider the equation 

L g $ = 4«5\ A + VV (2.13) 

This is treated as follows. Write 

$ = G + /i, (2.14) 

so that (2.13) gets replaced by 

L g h=l(l + G- 1 h)- 7 G- 7 p. (2.15) 
8 

In order to solve Equ. (2.15) we use the Green function G(x, x') with source at an arbitrary 
point x' G M. Define 

Q (x,x') =g ah {x'){x-x') a {x-x') b (2.16) 

Q(x, x') = Q {x, x') + ^g ab , c (x') (x - x') a (x - x'f{x - x'f + O°°(tf (x, x')) (2.17) 

in coordinate neighbourhoods of (x, x') E M x M, and extended to all of M x M as a 
smooth and positive function. Clearly we have 

G(x, x') = tt~ l l 2 (x, x') + ^ft + 0(Sll /2 (x, x')) (2.18) 

d x [G{x,x') - Q-^ix,^)] = 0(1), (2.19) 
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which relations also imply that 

G(x,x') = Q 1/2 (x-x') + 0(1) (2.20) 

d x [G(x, x') - Q 1/2 (x, x')} = 0(1). (2.21) 

It follows from the Appendix of Ref. [11] that there exists a unique positive solution h of 
(2.15), smooth on M \ A, bounded and with bounded first derivative on M. Thus 

h = 0(l), dh = 0(l) at A. (2.22) 

We now come to the construction of maximal initial-data sets (R 3 ,g, p) with R 3 arising 
as M \ A = S 3 \ A. Let p ab be a symmetric tensor, smooth on M \ A, which is trace- and 
divergence-free (TT-) with respect to g ab , where \i(g) > 0. Assume that 

Pab = O(fio 3/2 ), dp ab = 0(n- 2 ) at A. (2.23) 

Next solve Equ. (2.13), resp. (2.15) by taking p = p a bP ab ■ It then follows that, with 
9ab = ® 4 9ab, 

K[g\ = l -p ab r h on R 3 , (2.24) 
where p ab = §~ 2 p ab . Furthermore 

D a p\ = 0, p a a = 0. (2.25) 

Thus (R 3 ,^, p) is a maximal solution of the Einstein vacuum constraints. Setting 

m = /i+2h\ A (2.26) 

we have that (f := Qq 1 ^ 2 ) 
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in coordinates x a = VL~ l (x a + ^Tl c \\x b x c ). Thus (g,p) is asymptotically flat (in fact: 
asymptotically Schwarzschildian) with mass m and vanishing momentum. 

We now recall how TT-tensors p ab having the asymptotic behaviour (2.23) are con- 
structed. This is done on R 3 with metric g' ab given by 

g' ab = n~ 2 g ab , (2.30) 

which satisfies 

J* = ^ + (h) » ^ = °°° (k) ( 2 - 31 ) 



f?2 I ' vat> 

in the ^-coordinates. 

Again, by conformal invariance of the TT-condition we have to solve 

D'j>'\ = 0, p b = (2.32) 

with 

P'« b = O 1 (^) • (2.33) 
More specifically we want to find a TT-tensor p' ab with 

Pab = §(3M - r 2 5 ab ) + O 1 (^-) , (2.34) 



~c v"-™u-~o ■ ~auj i " i —3,, 

<yO \<y>0\0 

where C is a non-zero constant and x a = S ab x b . Note that (2.34) agrees asymptotically 
with (1.7), obtained for the SS metric. Using [14], this is accomplished as follows. Define 
the conformal Killing operator L, acting on covector fields A a by 

(L'X) ab = D' a \ b + D' b X a - 2 -g' ah {D' c \ c ) (2.35) 



and 



(D' ■ L'X) b := D' a (L'\) ab = A'X b + l -D' b {D' a \ a ) + K' b a \ b . (2.36) 



With our asymptotic conditions on g', the operator D' • L' is an isomorphism from the 
weighted Hdlder space C e fe+2 ' a (R 3 ) to C £ fc ; Q 2 (R 3 ) for all k G N and all < e < 1 (see [15] 

12 3 4 

for the precise definitions). We now seek fields Aa, Aa, Aa, Aa solving 

(D' ■ L'\) a = (2.37) 
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with 



L = /ia + 0°°f4) ( 2 - 38 ) 



L = W b + C°°(^) (2.39) 
L = e abc x b K c + 0°°(±-) (2.40) 



\ a = Qx a v + 0°°(lr) (2.41) 



where e abc is the Euclidean volume element, /i a , fi ab , k c and v are constants with fi ab 
symmetric and tracefree w.r. to 5 ab . This is done by observing that the leading terms in 
(2.38-41), say a a , have D' ■ L'a of fast decay, setting \ a = a a + 5\ a and solving 

D' ■ L'5X = —D' ■ L'a. (2.42) 

A 12 3 4 

Denote by Aa a basis of the 12-dimensional vector space spanned by Aa, Aa, Aa, Aa and 
define a tracefree tensor Q ab = Q( a b) by 

12 A 

Q ab = fJ2 c a(L' A) a6 , (2.43) 

A=l 

where / G C^°(R 3 ), non-negative and not identically zero. Then we seek constants ca 
solving the linear system 

J2E A BC B = -d A , (2.44) 

B 

where d,A is the vector (/i a = 0, fi ab = 0, K a = 0, v = YlixC) and 

Eab = E {A B) = I . f(L' \) ab (L' \) ab dV(g'). (2.45) 

It is shown in [11] that, when g' ab has no conformal isometry, Eab has trivial null space, 
and (2.44) has a unique solution. If g' ab is only conformally non-flat, E A b has at most a l- 

1 _ 2 

dimensional null space given by ca = (fJ> a , Hab = 0, k c , v — 0) where A a =A a (Hb)+ Aa (^ c ) 
is a conformal Killing vector. Clearly cIa is orthogonal to this null space so that (2.44) 
can still be solved, and solved uniquely if we require in addition that 

E caCa = 0. (2.46) 

A 
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Thus, given C in the definition of g^, we can find ca, whence Q ab - For this Q ab we solve 

(D' ■ L'W) a = D' b Q ab . (2.47) 
There is a unique solution W a to any equation of the form 

(D' ■ L'W) a = j a , (2.48) 
with j G C^ +2,a (R 3 ), lying in ,a (R 3 ). Writing the differential operator as it's flat-space 

o o 

analogue D • L plus the rest, we infer that 

(D -LW) a = ]a + 5j a = p a . (2.49) 

Since g' — 5 = (^^(l/f 2 ), the Euclidean components of p a are again in C^ +2,a (R 3 ). Now 
the arguments in [11] show that 

w x a P b x b + 7 'Par 2 e a bc x b L c 3x a M bc x b x c + QMa^r 2 Cx a , xu , fo 
W * = ~3 + ^ + ^ 2^ + ( 2 - 5 °) 

where 5W = O^^l/r 2 ^) with the constants involved in 5W = O^l/r 2 ^), d5W = 
0(l/f 3+e ), . . . d k 5W = 0(l/f 2+fc+£ ) depending only on pointwise weighted bounds on j 
and its first k derivatives. 

In the particular case where j a = D' b Q a b, and with the choices made by (2.45) and 
by solving (2.46,47), it was shown in [17] that P a , L ao and M ab are all zero and (7^0. 
Thus, by a computation, the tensor 

p' ab = Q ab - (L'W) ab , (2.51) 

in addition to being TT, satisfies (2.34). 



3 Critical Sequences 

We now consider critical sequences (CS's) g n of background metrics on M. We require that 
these metrics all have \\{g n ) > and tend smoothly to a metric g^ satisfying A 1 ((? OCJ ) = 0. 
This requirement of smooth convergence could be weakened considerably, but we shall not 
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attempt this. It is known that there are plenty of such sequences. In fact (see [16]), every 
metric g with \\(g) > can be deformed into a metric h with X\(h) < by changing it in 
an arbitrarily small subset of M, and then one could define g t — (1 — t)g +th, which is a 
continuous sequence of Riemannian metrics for t G [0, 1]. Since, by standard perturbation 
theory [17], Ai(<7 t ) depends continuously on t, there is a t G (0, 1) such that \i(gt ) = 
and Ai(<7t) > for t G [0,t ). Then define g n = (1 — t n )g + t n g to , where t n is any infinite 
sequence in [0, io) with lim^oo^ = to- This is a CS. 

It was shown in [10] that, along any CS g n , the time symmetric mass /i n in (2.7) goes 
to infinity. Furthermore there exists a constant E independent of n such that 

G n -n-^-^<^ n E^/\ (3.1) 

(We shall henceforth always denote positive, n-independent constants by E or E' with 
the understanding that they are not necessarily identical.) Also 

d(G n - 1T 1/2 ) < pL n E. (3.2) 

Note that Q, which because of Equ.'s (2.4,5) depends on n, also has an upper bound 
which is independent of n. Equ.'s (3.1,2) imply that 

G n - Q~ 1/2 < n n E (3.3) 

d[G n -Q- 1/2 ] < fi n E. (3.4) 

It also follows from [10] that 

G n - 1T 1/2 > fi n E, (3.5) 

whence 

G n - n^ 1/2 > fi n E. (3.6) 

All constants E, E' originate from the Schauder and Harnack inequality on open subsets 
of M and thus can be expressed in terms of global pointwise bounds on g n and a finite 
number of it's derivatives. It follows that these constants, in addition to being independent 
of n, can be taken to be independent of the choice of source point A. Thus, in the Green 
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functions G n (x,x'), the same bounds will hold with Q replaced by Q(x,x') and in turn 
by Q (x,x'), but with \i n now a function of x'\ /i n = /i n (x'). It is then easily shown [11] 
that, in fact, /i n (x') has a uniform bound from above and below in terms of /i n at any 
other point of M, say A. Writing now /i n for fJ> n \ A , we get 

E'Pq 1/2 (x, x') + < G n {x, x') < E[^ l/2 (x, x') + n n ] (3.7) 
d x [G n (x,x')-n 1/2 (x,x')} < Efi n . (3.8) 

Let N A be a fixed coordinate neighbourhood of A with chart y a centered at A. Let g s be 
a metric on M which coincides with the flat metric 5 ab in N A and let Vt l J 2 (x, x') be the 
Euclidean \y — y'\, extended smoothly as a positive function to all of M x M. Clearly we 
have 

ETl s (x, x') < Q {x, x') < EQ s {x, x'). (3.9) 

Thus 

E'[nj l/2 (x,x')+fi n ] < G n (x,x') < E[nj 1/2 (x : x')+pL n ] (3.10) 

d x G n (x,x') < D[nj\x,x')+fi n }. (3.11) 

Given the CS g n we can decompactify as shown in the last section, thus obtaining a 
sequence g' n of asymptotically flat metrics with components converging to a metric g'^ 
in the Holder norms C^ +3 ' a (R 3 ) for all k e N. For each g' n and for g'^ there is the 
operator D' ■ L', written down in Equ. (2.36) which maps C k+2,a (R 3 ) isomorphically onto 
C' ,a (R 3 ). Thus, if there is a sequence of sources j n converging to some joo in the C k,a - 
norm, the solution W n of (D' n ■ L' n )W n = j n converges to in the C k+2 ' a -iiorm. Hence 

A A 

the fields An, defined in the last section, converge to Aoo corresponding to the metric g'^. 
In particular the constants involved in the remainder terms in (2.38 - 2.41) can be taken 
to be independent of n. 

Let us now, for ease of presentation, suppose that the sequence g n and it's limit g^ 

A 

are "generic". By this we mean that no linear combination of the A's for these metrics 
is a conformal Killing vector. Then, by performing for each n the procedure described 
between Equ. (2.43) and Equ. (2.51), we obtain fields W n with P a = L a = M ab = and 
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fixed C ^ converging, in the +2,a -norms, to a with the same properties. Splitting 
the operators D' n ■ L' n into a flat-space part and the rest, it is standard to see [18] that 
the bounds on SW in (2.50) are uniform in n. When we now undo the decompactification 
we obtain, for each metric g n in the CS, a TT-tensor p ab on M obeying 

Pab = C 3Da{Ql ° /2)Db ^P - 9ablA + Q°°(Qo (3+£)/2 ), C>0 (3.12) 

with 0°° independent of n. We now have to solve, for each n, the equation (2.13) with 
P = PabP ab - This, in turn, requires to solve the integral form of equation (2.15), namely 

^^//^'fS* 1 ' (313) 

Suppose x G N\. Then the integral in (3.13) is dominated by the contribution coming 
from integration over x' G N A . Using (3.9, 10, 11, 12) and that 

E'Q 7 / 2 (x') < G~ 7 (x f ), (3.14) 

we see, that h is bounded by positive constants times the sum of the following integrals 

Mv) = [ , ( A + »nD) 7 ^rfy (3.15) 

J\y'\<R \y-y'\ V 12/ 1 / b l 

and 



B(y)= f i n f (A + -^dV, (3.16) 

•Vi<«\|j/| / Is/ 1 6 

where i? is a bound on the size of N A . It is elementary to check that 



A = 0(p- 3 ), B = 0{^) (3.17) 



with O independent of n. 

Similarly we find that dh is uniformly bounded by 

-7 



c (y) = [ , r^-r^ ivr\ + rm^ ( 3 - 18 ) 

J\y'\<R \y-y \ 2 \\y \ J \y \ b 



and 

-7 



D(y) = f i n [ (-L+f^D) r^d 3 y'. (3.19) 

J\y'\<R \\y I / |y l 
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We find that 

C(y)=0(^), D{y) = 0{^ 2 ). (3.20) 

Consequently 

h = Ofa 3 ), dh = 0{^ n 2 ). (3.21) 

The above estimates provide information on the quantities Q± associated with surfaces 
"f2 = small constant" as n gets large. Recall that 

e± = HTp ab n a n b . (3.22) 

The quantity H, in the present case, is given by 

O a O b \ ._ 1 

(3.23) 



H = D a n a = 



I Q l/2 $ ( ab _ ^\ Q + ^1/2$ flo 

4 ^ y n^v afe 



In A/a we can, for each n, change the y-coordinates centered at A to y a for which 5w>| A = 
5 ab and then change these, in N A \ A to (#,</?, r), where (0, </?, r) is related to y a like 
standard spherical coordinates to standard Euclidean coordinates. The quantity 

n a n b 9ab 



b = Att 

appearing in (3.23) can be expanded as 

b = 4 + hr + f 2 r 2 + f 3 r 3 + 0(r 4 ) (3.24) 

with fi, / 2 , fz smooth functions on (9, ip) e S 2 . 

Recall also that $ = G + h, G being the conformal factor giving the time-symmetric 
solutions conformal to the given background metrics g n . The quantity G, in turn, [19] 
can be written as a sum 

G = G loc + F, (3.25) 

where G loc is the Hadamard fundamental solution of L g in a neighbourhood of A, whence 
F is smooth and satisfies L g F = in this region. By [19] this region can be chosen 
independent of n. It is known [19] that G loc can be written as follows 

G loc = Q -i/2y + w ^ (326) 
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where V = 1 +0(Q), W\ A = and each of V and W permit an expansion in terms of s 2 , 
the (geodesic distance) 2 to the point A. Clearly s 2 can be expanded as s 2 = r 2 +gir z + . . . 
with g 1 , . . . smooth on S 2 , and the same is true for Q. Similar expansions hold for fl a and 
G loc . Making N A smaller, if necessary, inserting this into (3.23) and using (3.21) we find 
that the square bracket in (3.23) can be written as X + Y with 

X = 1 - + fi nai r 2 + fi n a 2 r 3 + 0(/i n r 4 ) + 0{^r) + 0(^ 2 r 2 ). (3.27) 

The coefficients a\ and a 2 are smooth on S 2 and remain bounded as n — > oo. We now set 

2 8ai 16a^ 64a? 

(3.28) 



8ai 16a 2 64a 2 

' n ^ 2 ^ 3 ~^ 3 

l^n f^n A*n 



and take ra large enough so that (#,</?, r n ) e A/a- Then Equ. (3.28) defines an embedded 
2-sphere S n C ATa- Evaluating X along E n , a computation gives that 



XI 



GO- 



(3.29) 



Thus the surface £„ becomes more and more minimal as n — > oo. The second term in 



Equ. (3.21) can be estimated from 



p ab n a n b 



4ft~ 1 / 2 $ 



-3 



(W) 1 / 2 

Now insert (3.12) into (3.30), using that 



(3.30) 



r = Ql /2 +O(Q ) = Ql /2 + 0(fi- 2 ), D a r = DA+O^ 1 ), 



We find that the square bracket in Equ. (3.30) is equal to 



r = i $- 3 [3 / i„c + o(^- £ )]. 



But, from (3.1), 
so that 

Thus, as n gets large, 
We have thus proved the 



$1 



3/i r 



+ 0(1), 



0+ < 0, 0_ > 



on S r 



D a rD a r = l+0( f i- 1 ). 

(3.31) 

(3.32) 
(3.33) 
(3.34) 
(3.35) 
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Theorem: Let g n be a generic initial sequence and (g n ,Pn) the unique initial-data set 
constructed for each n from the "free data" g n , p n with p n coming from the "same tensors" 
Q2b, so that p n is of the form (3.12) with a fixed constant C > 0. Then the 2-spheres S n 
given by (3.28) are future trapped for sufficiently large n. 

We call a CS special when, for each n and for g^, there is a linear combination of the 
A-fields defined in Sect. 2, which is a conformal Killing vector. (This has to be unique for 
goo and for g n with n large.) Then the above reasoning goes through unchanged, except 
that Equ. (2.44) is now solved with the side condition (2.46). Hence the above Theorem 
remains true with "generic CS" changed into "special CS" . 

The construction we have described for given g n and given function / appearing in 
(2.43), gives rise to a unique sequence of initial-data sets. There is of course much more 
freedom in the choice of extrinsic curvature. It is known [20] that, on the compact manifold 
M, there is an infinite-dimensional set of smooth TT-tensor r a b for any background 
metric g. Adding such a tensor to the one constructed from Q a b only changes p a b at 
higher multipole order and thus gives rise to a critical sequence for which the Theorem 
remains to be true. 

Finally, we could generalize our result by including angular momentum, that is to say 
construct data for which the W-field of Equ. (2.50) has no linear-momentum term P a 
and no quadrupole-contribution M ab , but both L a and C non-zero. The reason is that 
the angular momentum term in p a b gives, to leading order, a vanishing contribution to 
p a bn a n b ) so that p a bn a n b will asymptotically still have a sign for large n. 
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